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ABSTRACT 

The design of a robust compensator is considered 
for the SCOLE configuration using a frequency- 
response shaping technique based on the LQG/LTR 
algorithm. Results indicate that a tenth-order 
compensator can be used to meet stability- 
performance -robustness conditions for a 26th-order 
SCOLE model without destabilizing spillover 
effects. Since the SCOLE configuration is 
representative of many proposed spaceflight 
experiments, the results and design techniques 
employed potentially should be applicable to a wide 
range of large space structure control problems. 

Introduction 

Large space structures (LSS) have many properties 
that make them difficult to analyze and control 
[1]. They are mathematically modeled by computa- 
tionally difficult partial differential equations 
or high-order, lumped, ordinary differential 
equations obtained through finite element methods. 
LSS have many low and closely spaced resonant 
frequencies, a number of which typically fall 
within the controller bandwidth. In LSS, vibra- 
tional issues must be treated as a first-order 
effect; it is this characteristic of the LSS 
control problem that most distinguishes it from 
spacecraft control problems of the past. Addition- 
ally, inherent damping is low and/or improperly 
modeled. Coupled with stringent operational 
requirements for orientation, shape control, and 
vibration suppression, these properties present an 
unconventional and unresolved control design 
problem to the system analyst. 

A fundamental issue to be dealt with in any LSS 
control problem comes from the large amount of 
modeling error occurring in finite element models 
of such structures. In general, inaccuracy of 
modal data, such as elastic frequencies and mode 
shapes used to form coefficient matrices of the 
dynamic models, increases with increasing modal 
frequency. Hence, a frequency-dependent constraint 


is inherently imposed on the design process i\ that 
stabilization and performance requirements must oe 
met without allowing the input control energy to 
"spill-over" and excite and destabilize the lightly 
damped, poorly modeled high-frequency dynamics. 

At NASA's Langley Research Center, a LSS config- 
uration known as the Spacecraft Control Laboratory 
Experiment (SCOLE) was conceived for the purpose of 
evaluating and comparing large space structure 
control and identification concepts [2]. The SCOLE 
configuration (shown schematically In Figure 1) 
consists of a 130- foot flexible beam anchored at 
one end to the cargo bay of the space shuttle with 
an antenna reflector connected to the opposite end. 

The center of mass of the reflector is offset from 
the attachment point. The SCOLE configuration is 
representative of many proposed space flight exper- 
iments and space-based antenna systems. Control 
inputs are available from torque actuators located 
on the orbiter and force actuators at the reflector 
center. Attitude sensors are located at the 
reflector center. A typical SCOLE control task 
is to slew or change the line -of -sight of the 
antenna rapidly and damp any induced structural 
vibrations to the degree required for the precise 
pointing of the antenna. 

In this paper we consider a SCOLE large -angle 
slewing maneuver to have been completed and attack 
the problem of designing a model -based compensator 
to attenuate residual structural vibrational motion 
and antenna line -of -sight error. The SCOLE mathe- 
matical model is first discussed followed by des- 
criptions of the design objectives and the compen- 
sator design approach. Finally, results from the 
application of the design methodology to the SCOLE 
problem are presented. 

Mathematical Model 

The basic distributed-parameter mathematical model 
of the SCOLE configuration is described in [2], 
while nonlinear and linear ordinary differential 
equation models are found in [3] and [4], respec- 
tively. A linear finite-element model consisting 
of three rigid rotational modes and the first ten 
structural elastic modes is used in this study. A 
state -space realization of the modal model has the 
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form 


where 


and 


* F - Vf + b f u 

y F - Vf 


A f - diag( A^ A e ) 

A — [ ° ‘ 

\ 0 0 

6x6 


A e - diag( aJ. 4 A e ° ) 


A E“ 


-“i 


2x2 


( 1 ) 

( 2 ) 

(3) 

(<*) 

(5) 


( 6 ) 


for (i-1 , . . . , 10) . Equation (4) describes the rigid 
body contribution and equations (5) and (6) des- 
cribe the elastic contribution for ten vibrational 
modes of frequencies (i-1 10). A uni- 
form damping ratio of ■ f ■ 0.003, (i - 1 

10) , is assumed. The eigenvalues of A £ are given 
mathematically by 

and are shown in Table 1. 

Five control inputs are generated using three 
torque actuators (one per X,Y,Z axis) on the orbi- 
ter and two force actuators (X and Y directions in 
Figure 1) at the reflector center. Three attitude 
sensors (one per axis) are located at the reflector 
center. Sensor and actuator dynamics were not 
included in this study. Rigid-body inertias, 
mode shape and slope data from the finite element 
analysis combine to define the control effec- 
tiveness matrix Bp and output response matrix Cp. 


Design Objectives 


The design objectives of this study are to produce 
a multivariable, model -based, feedback compensator 
operating on attitude sensor data which will gener- 
ate force and torque Inputs to stabilize the rigid 
body modes; enhance the stability of lightly 
damped, low- frequency modes without destroying the 
stability of higher -frequency modes; meet pre- 
scribed closed- loop performance (bandwidth) speci- 
fications; and possess some degree of stability 
robustness to unmodeled dynamics. Since a low- 
order controller is sought, it was decided to 
employ full -state controller design with a reduced- 
order plant model. The full -order model is 
reserved for evaluation purposes. Order reduction 
for the design plant was performed using modal 
truncation. Past studies ([6], [7]) have indicated 
that a 0.1 rad/sec closed- loop performance band- 
width is sufficient to maintain antenna pointing 
control, and a design model composed of the rigid 
body plus the first three elastic modes in Table I 
is adequate to achieve this bandwidth. Higher 
bandwidths will typically require the addition of 
extra elastic modes to the design model. Denoting 
the 12th-order design model transfer matrix by 
Gp(s), for a unity-gain feedback compensator with 

transfer matrix G (s) , multivariable bandwidth will 
c 

be defined as the frequency below which the 
smallest singular value of the closed- loop response 
matrix 

G cl (s) - G p (s)G c (s)[ I + G p (s)G c (s) )' 1 (8) 

remains above unity for s - ju>. In our case, we 
seek a compensator such that 

o [ G __ ( j(*>) ] £ 1.0 for 0 £ w 5 0.1 (9) 

From the block diagonal structure of Ap in (1), the 
transfer matrix, Gp(s), of the 26th-order system 
may now be written as 


Analysis of (1) and (2) verifies that the system is 
completely controllable and observable. Attempts 
to reduce the number of control variables to the 
number of outputs retained controllability and 
observability but, in each three - control input 
combination, introduced 1 ightly- damped , low- 
frequency transmission zeros [5] into the model. 
Since the presence of such transmission zeros has 
been demonstrated to reduce system performance in 
large space structure controller designs [6], the 
compensator was designed with the original five 
inputs and three outputs. However, in order to 
avoid numerical ill-conditioning brought about by 
the different physical characteristics of forces 
and torques, the inputs were scaled so that the 

frequency response of the largest (cr(jw)) and 
smallest (a(jo>)) singular values of the transfer 

matrix of (1) and (2), denoted by Gp( jw) , were 

nearly equal at low frequencies (as shown in Figure 

2 ). 


G f (s) - G (s) + AG(s) (10) 

where AG(s) represents the transfer matrix of the 
remaining 14th-order (residual) modal system. In 
this form, the dynamics represented by AG(s) can be 
interpreted as an "additive perturbation" to the 
Gp(s) system and used as an approximate represen- 
tation of unmodeled dynamics for use in stability 
robustness tests. Specifically, it is established 
in [8] that the unmodeled dynamics AG(s) will not 
destroy the closed- loop stability so long as 

al G (s)[ I + G (s)G (s) J* 1 ) a[AC(s)] < 1 (11) 

c pc 

for all s - jo>, cj real. Condition (11) can be 
enforced in the compensator design stage to ensure 
that closed-loop stability will be preserved for at 
least that class of unmodeled dynamics whose 

spectral norm lies below a[AG(ju>)]. 
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The spillover effect on AG(s) due to the closed- 
loop compensation may be tested directly by 
applying a state-variable realization of G^Cs) to 

the full 26th-order model given by (1) and (2) and 
examining the eigenvalues of the composite system. 

A block diagram for the closed- loop configuration 
is shown in Figure 3. An approach for constructing 
G c (s) to stabilize G p (s) while satisfying (9) and 

(11) is presented in the next section. 


which, when used in an LQG fashion with the Kalman 
filter from Step 1 , asymptotically recovers the 
frequency response of the target loop gain over the 
low-frequency band. The resulting compensator is 
given by 

G (s) - F( sl-ft ) _1 H (17) 

c 

where 

F - B P (18) 

A T P + PA - PBB T P + qC T C - 0 (19) 


Compensator Design Approach 

The compensator design approach employed to meet 
the foregoing design objectives can be viewed as a 
variation of the well-known Linear -Quadratic - 
Gaussian/Loop-Transfer-Recovery (LQG/LTR) algorithm 
([9), [10]). In the standard LQG/LTR approach, 
with the loop in Figure 3 broken at the output, a 
Kalman filter (G K p) is designed to meet the 

complete set of stability-performance-robustness 
objectives. Thereafter, an optimal linear regu- 
lator is constructed such that the composite LQG 
compensator (G^) loop gain behavior asymptotically 

approaches (recovers) that of G^ in the sense that 


k - A-BF-HC (20) 

Step 3 

Attempt to adjust q in Step 2 until the desired 
bandwidth condition (condition (9)) is met. Also 
check stability robustness by (11). If an exces- 
sively high q (indicated by violation of (11)) is 
required to achieve the required bandwidth, turn 
down the Kalman filter gain (by increasing p in 
(14)) to "loosen" the target loop. In effect, this 
procedure reduces the target bandwidth until satis- 
faction of (11) is possible. The final design is 
accomplished by Iteratively adjusting the linear 
regulator and Kalman filter design parameters until 
an appropriate compromise is made between bandwidth 
and stability robustness. 


G p <j, )G c ( ju>) — G KF (ju>) 

pointwise in w. Direct application of this LQG/LTR 
procedure to large space structures problems 
results in extremely conservative designs which 
cannot meet reasonable performance specifications 
[6]. However, the LQG/LTR structure still provides 
a viable approach for model-based controller 
synthesis when the standard procedure is modified 
in the following manner. 


Step 1 

Denote a state-variable realization of G p (s) by 

x - Ax + Bu (12) 

y - Cx (13) 

Select the design parameters L and p in the Kalman 
filter algorithm 

AQ + QA T + LL T - — QC T CQ - 0 (14) 

H - — QC T (15) 

such that 

g kf ( s) “ c( SI ‘ A >’ 1h (16) 

achieves a desired (target) loop gain for 
G (s)G (s) over some low-frequency band containing 
P c 

the design bandwidth. 

Step 2 

By successively increasing q > 0 in equation (19) 
(to follow), design an optimal linear regulator 


In large space structures applications, the 
inability (at Step I) to meet loop gain magnitude 
over the desired bandwidth or (In Step 3) the 
production of too small a compromise bandwidth can 
often be overcome by the inclusion of additional 
flexible modes Into the design model [6]. 

If, as In the SCOLE application to follow, an 
order- reduction study is performed on the resulting 
compensator, the complete set of stability- 
performance -robustness conditions needs to be 
re-evaluated with the reduced- order compensator. 


SCOLE Application 


Figure 4 shows the frequency response of the 12 th- 
order (LQG) compensator, G c , resulting from an 


application of the foregoing procedure to the 12 th- 
order SCOLE design model, G . The figure indicates 


a well-behaved lead- lag structure with a 20 
db/decade roll-off. Eigenvalues of the corres- 
ponding A-HC, A-BF, and A matrices are given in 

Table II. The frequency response of G in equa- 

CL 


tion (8) with the 26th-order evaluation model used 
in place of G is shown in Figure 5. Figure 5 

demonstrates the satisfaction of the 0.1 rad/sec 
bandwidth requirement and a 60 db/decade roll-off. 
Eigenvalues of the composite system resulting from 
the LQG compensator applied to the evaluation model 
are given in the first column of Table III. The 
data Indicate that the compensator was designed to 
concentrate on stabilizing the rigid body modes 
(first three entries of the A-HC block) without 
disturbing the stability of the three elastic modes 
(last three entries of the A-HC block) of the 
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design model. Table III also shows that there is 
insignificant spillover into AG. The stability 
robustness test (11) shown In Figure 6 shows more 
than 10 db robustness margin. 

The possibility of a reduced-order compensator 
(ROC) satisfying the design conditions was also 
investigated. The methods of balanced realization 
[11], Hankel-norm reduction [12] and a method based 
on stable factorization [13] were employed. The 
ratio of largest to smallest Hankel singular values 
was 0.016 so little reduction based on nearly 
uncontrollable or unobservable compensator modes 
was expected. All of the methods gave similar 
results. In each order reduction method, only a 
lOth-order ROC would stabilize the design model. 

The stable factorization results were judged to be 
(slightly) better and will be discussed herein. 

Table II shows the eigenvalues o£ the corresponding 
A matrix in the ROC, denoted by A^^ . Figure 7 

shows a frequency response of the ROC. A compar- 
ison of Figures 4 and 7 shows that the only 
difference between the LQG and ROC frequency 
response plots is the removal of the dip in a at 

the frequency of the third elastic mode. The 
importance of this characteristic can be seen from 
Figure 8 where the robustness condition (11) is 
evaluated using the ROC in place of G^ . Figure 8 

indicates that an effect of the-order reduction is 
a reduction of stability margin at the frequency of 
the third elastic mode. A more positive effect 
from the reduced order compensation is seen in the 
second column of Table III where the eigenvalues of 
the ROC applied to the evaluation model are presen- 
ted. The stability of the rigid-body modes from 
the LQG compensator is preserved with the auxiliary 
effect of adding stability to the first elastic 
mode. The net effect of the ROC is to enhance the 
stabilizing effect of the LQG compensator at the 
expense of a reduction of stability robustness 
margin. 

Concluding Remarks 

A loop- shaping procedure similar to that used in 
the LQG/LTR approach was used to design a model - 
based compensator for the SCOLE configuration, a 
generic large space structure configuration 
conceived for the purpose of evaluating and 
comparing control and identification approaches. 
Initially, the inputs of a full 26th-order SCOLE 
model were scaled to avoid numerical difficulties. 

A 12th-order controller design model was afterwards 
constructed from the full -order model using modal 
truncation. Applying a modification of the LQG/LTR 
technique to the design model produced a 12th- order 
model-based compensator satisfying stability- 
performance -robustness design conditions. Finally, 
an order- reduction technique based on stable 
factorization was used to produce a lOth-order 
compensator for controlling the full-order model 
without destabilizing spillover effects. It was 
noted that order reduction can have beneficial 
effects on closed- loop stability but may reduce 


stability robustness margins. Since the SCOLE 
configuration is representative of many proposed 
spaceflight experiments, the results and design 
techniques employed should potentially be 
applicable to a wide range of large space structure 
control problems. 
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TABLE I 

EIGENVALUES OF A_ 
E 


Mode 

Eigenvalue* 

1 

(-0.00524, 1.747) 

2 

(-0.00591, 1.970) 

3 

(-0.0513, 5.108) 

4 

( -0 .0224 , 7.449) 

5 

(-0.0387, 12.903) 

6 

(-0.0898, 29.925) 

7 

(-0.104, 34.657) 

8 

(-0.232, 77.165) 

9 

(-0.243, 80.993) 

10 

(-0.446, 148.780) 


* TABLE II 

EIGENVALUE ANALYSIS OF COMPENSATORS 


A-HC 

A-BF 

k 


^SlOC 

(-0.126, 0.126) 

(-0.0563, 

0.101) 

(-0.420, 

5.125) 

| (-0.776, 

1.876) 

(-0.126, 0.126) 

(-0.0834, 

0.0837) 

(-1.038, 

2.027) 

(-0.126, 0.126) 

(-0.0959, 

0.104) 

(-0.203, 

1.974) 

(-0.204, 

1.973) 

(-0.00524, 1.747) 

(-0.178, 

1.976) 

(-0.0489, 

0.184) 

(-0.0500, 

0.185) 

(-0.00591, 1.970) 

(-0.905, 

1.956) 

(-0.197, 

0.237) 

(-0.197, 

0.237) 

(-0.0153, 5.108) 

(-0.418, 

5.130) 

(-0.209, 

0.210) 

(-0.209, 

0.210) 


^ TABLE III 

EIGENVALUES OF FULL- ORDER SYSTEM FORCED 
BY LQG AND REDUCED -ORDER COMPENSATORS 


LQG 

ROC 

(-0.0569, 0.101) 

(-0.0550, 

0.0998) 

(-0.0834, 0.0837) 

(-0.0834, 

0.0837) 

(-0.0971, 0.103) 

(-0.0960, 

0.103) 

(-0.179, 1.976) 

(-0.181, 

1.975) 

(-0.905, 1.975) 
(-0.418, 5.130) 

| (-0.634. 

1.823) 

(-0.122, 0.128) 

(-0.122, 

0.129) 

(-0.126, 0.126) 

(-0.126, 

0.126) 

(-0.126, 0.125) 

(-0.128, 

0.129) 

(-0.00524, 1.747) 

(-0.152, 

1.752) 

(-0.00591, 1.970) 

(-0.00603 

, 1.970) 

(-0.0153, 5.108) 

(-0.0137, 

5.109) 

(-0.0243, 7.449) 

(-0.0242, 

7.449) 

(-0.0383, 12.903) 

(-0.0383, 

12.903) 

(-0.0898, 29.926) 

(-0.0898, 

29.259) 

(-0.104, 34.657) 

(-0.104, 

34.657) 

(-0.232, 77.165) 

(-0.231, 

77.165) 

(-0.243, 80.993) 

(-0.243, 

80.992) 

(-0.446, 148.780) 

(-0.446, 

148.780) 


Eigenvalues presented in 
(Real, + Imaginary) format. 
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